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' Abstract. The yields of charmed hadrons emitted by strangeness rich QGP are 

^Sj ' evaluated within chemical non-equilibrium statistical hadronization model, conserving 

strangeness, charm, and entropy yields at hadronization. 



1. Introduction 

I A relatively large number of hadrons containing charm and bottom quarks are expected 



to be produced in heavy ion (AA) collisions at the Large Hadrons Collider (LHC). 



O 

cn 

lO ! Because of their large mass c, c quarks are produced predominantly in primary parton- 

^ I parton collisions pp, at RHIC and thus even more so at LHC. These heavy flavor 

O I quarks participate in the evolution of the dense QCD matter from the beginning. In 

' view of the recent RHIC results it can be hoped that their momentum distribution could 

P-i' m 

I I reach approximate thermalization within the dense QGP phase j2]. 

We study in this report heavy quark hadronization, and we are particularly 
^ ' interested in how strangeness influences heavy and multi-heavy hadron yields 
(containing more than one heavy quark). In the presence of deconfined QGP multi-heavy 
^ I hadrons are formed from heavy quarks created in the initial NN collisions. Therefore 
ci I yields of these hadrons are expected to be enhanced as compared to yields seen in single 
NN collisions jHEj- This type of enhancement can be considered to be an indicator of 
the presence of deconfined QGP phase, for reasons which are analogous to those given 
in the context of multi-strange (anti) baryon enhancement jH]. J/'^{cc) yields were for 
this reason obtained in the kinetic formation and dissociation model |31 Hj but without 
consideration of a strangeness content in the QGP phase. 

Differing from other recent studies which assume that the hadron yields after 
hadronization are in chemical equilibrium |H], we obtain the charm hadron yields 
in the statistical hadronization approach based on an abundance of m, rf, s quark pairs 
fixed by the bulk properties of a chemically equilibrated QGP phase. Since we do 
not yet know to sufficient precision the heavy flavor production at LHC or even RHIC 
energies, we present particle yields normalized by the total charm flavor yield, such that 
the particle yield results we consider are little dependent on the unknown total yield of 
charm. The results we present were obtained for dNc/dy = c = 10. 
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2. Fast Hadronization 

We study within the statistical hadronization model (SHM) the charmed hadron 
production at fixed strangeness pair yield per entropy s/S. This in general requires 
the chemical non-equilibrium version of the SHM. The important parameters of the 
SHM, which control the relative yields of particles, are the particle specific fugacity 
factor A and space occupancy factor 7. The fugacity is related to chemical potential 
fi = Tin A. We take Aj = 1 for all fiavors since the small particle-antiparticle yield 
asymmetry at RHIC and LHC is not relevant here. 

The ratio of the number of given type particles produced, to the number expected 
in chemical equilibrium is (nearly) the phase space occupancy 7, which is the same for 
particles and antiparticles of the type considered. We use occupancy factors 7,^ and 
7^ for QGP and hadronic gas phase respectively, tracking every quark fiavor {i = q, 
s, c). For each hadron, the phase space occupancy 7^ is the product of 7^'s for each 
constituent quark. For example for the charmed meson D(cg), 7)^ = 7^7^. The value 
of 7^ has upper limit arising from the condition of Bose- Einstein condensation of pions 

The number of particles of type 'i' with mass rui per unit of rapidity is: 

^ = 7.nr— . (1) 
dy * dy 

Here dV/ dy is system volume associated with the unit of rapidity, and ri^^ is a Boltzmann 
particle density in chemical equilibrium: 

nr = 9^1 ^A,exp(-VPM^/T) = \^^g,{m,/TYK,{m,/T). (2) 

During a fast transition between QGP and HG phases strange and heavier quark 
fiavor yields are preserved, as are the entropy per unit of rapidity [Hj, and the specific, 
per rapidity, hadronization volume: 

dNl_ _ dNl_ _dNi ■_ ^ ^. dS^ _ dS^ _ dS__ dV^ _ dV^ 
dy dy dy ' ' ' dy dy dy ' dy dy 
The yields of hadrons after hadronization are given by Eq. (^, the three unknown 
7^, 7^ and 7^ can be determined from their values in the QGP phase, or equivalently 
corresponding fiavor yields, i.e. 7^^ or dN^/dy given in Eq. 0. 7^ 7^ 1 implies that 
hadron yields are in general not in a chemical equilibrium. We will show how this 
infiuences the relative yields of heavy fiavored particles in the final state. 
The number of strange quark pairs determines the value of 7^: 

dNg dV r H / H eq , H2 eq\ , H2/o H cq , eq n , q H3 eq] /.n 

~dy^ly^^' (7^ + % «y j + 7. [2% + ^shid) + 37. , (4) 

where n'^'^ are densities of strange mesons and baryons calculated using Eq. (Q) in 
chemical equilibrium, n^^id = ''^^^ + Psn"^^ and Pg is the strangeness content of the 
T). The pattern of this calculation follows an established approach by using SHARE 1.2 
program JU] in the calculation. 
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For charm hadrons is obtained from: 
dNr dV 



dy dy 



7?<. + if ' «lid + K<% + 27?n^?J ; (5) 



op ' / c V c hid ' I q ccq 



where open 'op' charm yield is: 

<p = ifn^ + ifnL + 7? ' + ifl^n^^, + 7? ' n^,; (6) 

Here and densities of D and Dg mesons respectively in chemical equilibrium, 

n^^g is equilibrium density of baryons with one charm and two light quarks, ra^^^ is density 
of baryons with one light, one s and one c quarks, n^^^ is density of baryons with one 
charm and two strange quarks (f2°(ssc)) and ^chid is the density of particles with both 
a charm quark and an anticharm quark (C=0, S=0). 

The use of the hadron phase space (denoted by H above) does not imply the presence 
of a real physical 'hadron matter' phase: the SHM particle yields will be attained solely 
on the basis of availability of this phase space. In fast hadronization one can assume 
that there are practically only free-streaming particles in the final state. 

Thinking in these terms, one can imagine that especially for heavy quark hadrons 
some particles are pre-formed in the deconfined plasma, and thus the heavy hadron 
yields may be based on a value of temperature which is higher than the global value 
expected for other hadrons. For this reason we will study in this work a rather wide 
range 140 < T < 260 MeV. 

We will use as a reference a QGP state with the ratio strangeness (s) to entropy (S) 
fixed. We take it be in the range s/ S =0.03 - 0.04, and consider dV/dy =600 - 800fm^ 
at T = 200 MeV. This corresponds to total particles' multiplicity of about 4000-5000 
after hadronization for T=140 MeV, according to calculations using SHARE 1.2 |10j . 
This is within the expected range of total hadron multiplicity per unit of rapidity for 
LHC We compare our results to the 'benchmark' reference yield which is obtained using 
7? = 7^^ = 1 (chemical equilibrium). 



3. Yields of charmed hadrons 



3.1. D, Dg meson yields 

Considering Eq. ((H), and using 7^, 7^^, 7^ at a given T we have now all the inputs 
required to compute relative particle yields of all charmed hadrons. In figure ^ we 
show in the upper panel the fractional yields of non-strange D/Nc and strange Dg/Nc 
charmed mesons normalized by the total number of charm quarks as functions of 
hadronization temperature. The dashed lines are for chemical equilibrium. The extreme 
upper and lower lines are for s/S = 0.03 (purple, solid and dash-dot lines), while the 
central lines are for s/S = 0.04 (green, solid and dash-dot lines). We recall that the 
chemical equilibrium of strangeness in QGP corresponds to s/S ~ 0.03. The yield of 
D is the sum over all (qc) states i.e. D° and {D~^ ^ D^) mesons and resonances. 
Similarly, Ds{sc) is the sum over all resonances of Dg mesons. 
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Figure 1. The upper panel: charm meson yield per total charm yield A'^c- 
Equilibrium (dashed lines, blue ), non-equilibrium (solid and dash-dot lines are 
for Dg and D, respectively). The lines for s/S = 0.04 (green) are more centrally 
located in each frame, outer (purple) lines are for s/S = 0.03. The lower panel: 
ratios D/Dg as a function of T. Solid line is for s/S = 0.04, dash-dot line is for 
s/S = 0.03, dashed line is for 7^ = 7^^ = 1 (chemical equilibrium). 

In the lower panel we show the ratios D/Ds as a function of hadronization 
temperature. There is a considerable deviation in the ratios obtained for a fixed s/S 
from the chemical equilibrium results, except for the accidental values of T where 
the equilibrium results (dashed line) cross the fixed s/S results. While the chemical 
equilibrium results show significant difference between strange and non-strange heavy 
mesons, for the fixed s/S case, with increasing s/S and towards low T these yields 
become nearly equal. If the heavy meson yields are established at temperatures similar 
to regular hadrons, relative enhancement by factor 2-2.5 in strange-heavy mesons is to 
be expected for LHC. For RHIC conditions {s/S = 0.03, c = c = 2) D/Dg ratio is lower 
than the chemical equilibrium value for T < 0.2. This ratio decreases by factor of 1.6 
for T ^ 140 MeV. 
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Figure 2. cc/N^ relative yields as a function of hadronization temperature T 
for s/S = 0.03 with reference volume dV/dy = 600 fm^^ at T = 200 MeV 
(upper panel); and s/S = 0.04 with dV/dy = 800 fm-^ at T = 200 MeV(lower 
panel). Dashed lines are for chemical equilibrium. 



3.2. Yields of hadrons with two heavy quarks 

The hadron yields with two heavy quarks are model dependent, being proportional to 
1/dV/dy because 7/^ for heavy quarks is proportional to 1/dV/dy. Thus our result is 
dependent on the reaction volume, or on the total assumed charm yields. In figure |21 
we show the yield of hidden charm cc (sum over all states of cc) mesons normalized 
by the square of charm multiplicity A*"^ as a function of hadronization temperature T 
for two different reference values of volume at T = 200 MeV: dV/dy = 600 fm^ (upper 
panel), and dV/dy = 800 fm^ (lower panel). We consider again cases with s/S = 0.03 
(upper panel, solid line) and s/S = 0.04 (lower panel, solid line). For comparison, the 
chemical equilibrium cc mesons yields are shown (dashed lines on both panels). Entropy 
per unit of rapidity dS/dy is conserved, during cooling/expansion of the QGP, which 
imphes VT^ ~Const. In the second case considered rapidity density of strangeness is 
78% higher than in the first. 

The yield of cc mesons is much smaller at fixed s/S = 0.04 than in chemical 
equilibrium for the same dV/dy for large range of hadronization temperatures. For 
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s/S = 0.03 the effect is similar, but charmonium suppression is slightly less pronounced. 
This new mechanism of charmonium suppression occurs due to competition with the 
yield of strange-heavy mesons. The enhanced yield of Ds in effect depletes the pool 
of available charmed quark pairs, and fewer hidden charm cc mesons are formed. For 
particles with two heavy quarks the effect is larger than for hadrons with one heavy 
quark and light quark(s). 

4. Conclusions 

We have considered here the abundances of heavy flavor hadrons within the statistical 
hadronization model. While we compare the yields to the expectations based on 
chemical equilibrium yields of light and strange quark pairs, we present results based 
on the hypothesis that the QGP entropy and QGP flavor yields determine the values of 
phase space occupancy 'yf,i = q, s, c, which are of direct interest in study of the heavy 
hadron yields. 

We studied how the (relative) yields of strange and non-strange charmed mesons 
vary with strangeness content. For a chemically equilibrated QGP source, there is 
considerable shift of the yield from non-strange D to the strange Dg. Since the expected 
fractional yield Dg/Nc — 0.2 when one assumes 7^ = 7^ = 1, one should be able to 
falsify chemical (non-)equilibrium hypothesis easily: the expected enhancement of the 
strange heavy mesons being at the level of 40% when s/S = 0.03, and yet greater, when 
greater strangeness yield is available. 

Another consequence of this result is that we find a relative suppression of the multi- 
heavy hadrons, except when they contain strangeness. The somewhat ironic situation 
is that while higher charm QGP yield enhances production of cc states, the enhanced 
strangeness suppresses this effect. 
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